Phase Transitions in Trajectories of a Superconducting Single-Electron Transistor 

Coupled to a Resonator 



<N 

o 

<N 



O 

CO 

<3 

> 
<N 

<N 

• 

O 



Sam Genway, Juan P. Garrahan, Igor Lesanovsky, and Andrew D. Armour 
School of Physics and Astronomy, The University of Nottingham, Nottingham, NG7 2RD, United Kingdom 

(Dated: May 31, 2012) 

Recent progress in the study of dynamical phase transitions has been made with a large-deviation 
approach to study trajectories of stochastic jumps using a thermodynamic formalism. We study 
this method applied to an open quantum system consisting of a superconducting single-electron 
transistor, near the Josephson quasiparticle resonance, coupled to a resonator. We find that the 
dynamical behavior shown in rare trajectories can be rich even when the mean dynamical activity 
is small and thus the formalism gives insights into the form of fluctuations. The structure of the 
dynamical phase diagram found from the quantum-jump trajectories of the resonator is studied 
and we see that sharp transitions in the dynamical activity may be related to the appearance and 
disappearance of bist abilities in the state of the resonator as system parameters are changed. We also 
demonstrate that for a fast resonator, the trajectories of quasiparticles are similar to the resonator 
trajectories. 

PACS numbers: 05.70.Ln, 03.65.Yz, 42.50.Lc, 85.85.+j, 85.35.Gv 



I. INTRODUCTION 

Transitions in the dynamical behavior of open quan- 
tum systems may be seen for a variety of driven systems. 
Notable examples include the laser [l|, atoms coupled 
to optical cavity modes such as in the micromaser 0, Q 
and pumped Bose-Einstein condensates mechanical 
resonators coupled to optical cavities [5|, and nanoelec- 
tromechanical systems [6]. Recently, an alternative per- 
spective on the dynamics of open quantum systems has 
been developed [7|, [8[ by studying ensembles of quantum- 
jump trajectories. Applying the so-called s-ensemble 0- 
[11| has given insights into dynamical crossovers by ex- 
tending the space of parameters @, H, EH . The new 
parameter V is a counting field which may be adjusted 
to give ensembles of trajectories biased towards fewer or 
more quantum jumps. Varying s allows the phase struc- 
ture of quantum-jump trajectories to be explored in a 
way analogous to tuning parameters across an equilib- 
rium phase transition, such as the temperature in a fer- 
romagnetic transition. 

The s-ensemble is developed in a formalism using a 
large-deviation method [14J, which treats the statistics 
of dynamical trajectories in the same way equilibrium 
statistical mechanics treats statistical ensembles of con- 
figurations [9I-HH. Just as order parameters, such as the 
magnetization for a ferromagnetic transition, are used 
to characterize configurations, we can use a dynamical 
order parameter to characterize quantum-jump trajecto- 
ries. The number of quantum jumps K in a given time t 
is the extensive conjugate variable to the counting field s. 
We therefore consider the activity fc, defined as the num- 
ber of quantum jumps in unit time K/t, as a dynam- 
ical order parameter. Characterizing dynamical phases 
in this way has been shown to give rich phase diagrams 
with phase boundaries marked by discontinuities in fc, 
or its derivative with respect to system parameters, oc- 
curring at first and second-order transitions [8] . For the 



unbiased physical dynamics when 5 = 0, the transitions 
typically appear rounded and we understand these to be 
crossovers in dynamical behavior, which become transi- 
tions in an appropriate thermodynamic limit [l5[ . 

Physical dynamics correspond to the s = ensemble 
of trajectories. The method gives a new perspective on 
these unbiased s = dynamics through the structure 
of phase diagrams parametrized by s. Furthermore, it 
should be noted that the statistics of trajectories when 
s / can be inferred from accurate measurements of the 
physical trajectories. The s-ensemble has already shown 
interesting results for driven 2- and 3-level dissipative 
systems. For example, it has been shown that for a par- 
ticular ratio of driving and dissipation Q, the two- level 
system exhibits a surprising self-similarity in the fluc- 
tuation properties for all s, so that all trajectories look 
the same if rescaled by the number of counts, k. For 
the three-level system studied in [7] it was shown that 
the intermittent bursts of quantum jumps [16| in the dy- 
namics arise from a sharp crossover in activity occurring 
at s = 0. In Ref. [8], dynamical phases of the micro- 
maser were studied. The micromaser consists of a series 
of two-level atoms passing through a microwave cavity 
in their excited state. As the flux of atoms is increased, 
the atoms drive the cavity through different dynamical 
states. It was shown that, as a function of the counting 
field s and the atom flux, distinct dynamical phases are 
seen, characterized by different values of the activity. 

In this paper, we examine a more complex system con- 
sisting of two coupled components with non-trivial inter- 
nal dynamics. We study a system constructed by cou- 
pling a resonator to a superconducting single-electron 
transistor (SSET) close to a Josephson quasiparticle res- 
onance (JQP) [17H24]. Both the resonator and SSET 
are themselves open quantum systems because the res- 
onator loses energy to its surroundings and Cooper pairs 
on the SSET break into quasiparticles which, one-by- 
one, decay incoherently off the SSET into a second lead. 
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This allows us alternative ways of generating trajecto- 
ries; we may either keep a time record of quanta emitted 
by the resonator or measure charges leaving the SSET 
island. While SSET-resonator systems were first demon- 
strated with nanomechanical beam resonators [23], su- 
perconducting stripline resonators have also been used in 
more recent experiments [24, 29]. It is the latter set up 
which, along with very recent progress towards detec ting 
single microwave photons from stripline resonators |25l - 
[28[, motivates our study of quantum jump trajectories 
presented here. 

This paper is organized as follows. In section [III we 
introduce the open-system model for the SSET and res- 
onator before describing how the trajectories of the sys- 
tem may be studied using the s-ensemble. In section HTH 
we present numerical results for the trajectories of quan- 
tum jumps in the resonator. Then, in section HV1 we dis- 
cuss a mean-field theory which describes the resonator 
limit cycles at s = 0. From this we develop an effective 
stochastic master equation for the resonator, from which 
we derive the activity of the resonator for all values of s. 
Finally, in section[V]we demonstrate briefly that studying 
the trajectories of quasiparticles allows the resonator tra- 
jectories to be inferred in certain regimes. In section IVT1 
we state our conclusions. 



II. MODEL AND FORMALISM 

The tandem of a resonator weakly coupled to an SSET 
near the JQP resonance may be explored with an open 
quantum system description [6j. The SSET island con- 
sists of a left- and a right-hand Josephson junction, a 
gate capacitor across which a gate voltage is applied 
and a capacitively-coupled resonator, shown in Fig.[TJa). 
The detuning from the JQP resonance in the left-hand 
Josephson junction can be controlled with drain-source 
and gate voltages, Va s and V g , while a resonator gate 
voltage also allows the strength of the coupling to the 
resonator to be adjusted. Near the JQP resonance, 
Cooper pairs tunnel coherently across a Josephson junc- 
tion, causing oscillations between SSET island states |0) 
and 1 2), with zero and two extra charges on the SSET 
respectively. However, the incoherent tunneling of quasi- 
particles through the right-hand lead takes the SSET is- 
land to the state |1), with a single extra charge, and 
then back to the state |0), as sketched in Fig. QJb). The 
coupling of a resonator to the SSET island allows for en- 
ergy exchange between resonator and SSET. Depending 
on the sign of the detuning from the JQP resonance, the 
resonator is either driven into states of self-sustaining os- 
cillation , or experiences a cooling effect [23| because 
of its coupling to the SSET. 

In the open-system model, the coherent oscillations of 
Cooper pairs across the left-hand Josephson junction, the 
motion of the resonator and the coupling of the resonator 
position to the charge on the SSET island are all included 
in an effective Hamiltonian. While the resonator may be 
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FIG. 1. (a) Schematic diagram of the SSET-resonator sys- 
tem showing left and right Josephson junctions labeled 'L' 
and 'R', with gate voltage V g applied across a gate capacitor 
and drain-source voltage Vds across both Josephson junctions. 
The resonator is also capacitively coupled to the SSET island, 
(b) Illustration of the JQP cycle. Cooper-pair tunneling oc- 
curs at the left-hand junction between states |0) and |2). This 
is interrupted by the quasiparticle decay to state |1) which in 
turn decays back to state |0). 



a stripline resonator [24( or a nanomechanical beam [23| , 
for clarity we will stick to language more appropriate 
for a mechanical oscillator. Both the dissipative effects 
of the environment on the resonator and the decay of 
quasiparticles from the SSET island and the right-hand 
lead can be described by Lindblad operators such that the 
dynamics of the SSET-resonator system may be found 
from a master equation 0, [30| . 



A. Master Equation 

We will focus on the regime where the SSET island 
exhibits coherent oscillations of Cooper pairs across the 
left-hand junction between charge states of the SSET is- 
land |0) and 1 2). Including the quasiparticle decay to 
the right-hand lead requires the inclusion of a third is- 
land state |1) with just a single quasiparticle. Treat- 
ing the resonator as a single-mode harmonic oscillator 
with frequency oo and considering the fluctuations in the 
resonator displacement to be small compared with the 
SSET-resonator distance, allows the coupling to be well 
approximated by expansion to linear order in the oscilla- 
tor displacement. The full Hamiltonian for the coherent 
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dynamics is then [6, 30]. 

H =A£|2)(2| - ft (|0)<2| + |2)(0|) + Wo 

+ Ca; i (at+o)(|l>(l| + 2|2)<2|). (1) 

The difference in electrostatic energy between states |0) 
and 1 2) is given by AE and Ej is the Josephson energy of 
the junctions. The displacement x s is the change in the 
equilibrium position of the resonator when one charge 
is added to the SSET island, with C = (Huj 3 m/2) 1 ^ 2 . 
In addition to the coherent dynamics governed by the 
Hamiltonian, dissipation occurs due to both quasiparticle 
tunneling off the SSET island and the coupling of the 
resonator to its environment. Including these effects, the 
full master equation reads 

p = W(p) = ~[H,p] + £ q . p . (p) + C osc . (p) (2) 

where p is the density matrix for the combined SSET and 
resonator system. The quasiparticle decay from states |2) 
to 1 1) and from |1) to |0) is described by 

£ q . p .(p) =r(|l)(2|p|2)(l| + |0)(l|p|l)(0|) 

-E({|2)(2|,p} + {|l)(l|,p}) (3) 

where {•, •} denotes the anticommutator and T is the 
quasiparticle decay rate which is, for simplicity, taken to 
be the same for the two processes |2) — » |1) and |1) — >• |0). 
The environment of the resonator will also be included 
via 

Cosc. (p) = 7ext«P« f - {a f a, p} (4) 

where 7 ext is the rate, per energy quantum Hlj, at which 
the resonator loses energy to its environment. The Lind- 
blad form (J2J) describes a thermal bath at zero tempera- 
ture which, although hard to achieve with a mechanical 
resonator, is valid for the higher- frequency stripline res- 
onators where fc#T <C fux). Indeed, Eq. (J2j) is valid for a 
broad range of oscillator frequencies, provided the cou- 
pling between the SSET and the resonator is weak. (A 
summary of the practical requirements for a system to be 
well-described by this model is given in the Appendix.) 

We assume throughout the model parameters [l8| T = 
Vds/eRj, with Rj = h/e 2 the SSET junction resistance. 
We can define the coupling strength in terms of the di- 
mensionless parameter A, where A 2 = mu 2 x 2 /eVa s will 
be much smaller than unity for the weakly-coupled sys- 
tem studied in this paper [3l|. 

A characteristic of this system is the ability of the 
SSET to drive @, or cool 0, H, [H the resonator 
depending on the sign of AE, which is set by the de- 
tuning of the SSET bias point from the JQP resonance. 
When AE < 0, the SSET can drive the resonator into 
different steady states as the SSET-resonator coupling 
strength is increased. When decoupled from the SSET, 
the resonator remains in its ground state. Upon increas- 
ing the coupling, a continuous crossover is reached where 
the oscillator enters a state of self-sustaining oscillation. 



This is illustrated by the distribution over the number 
states P(n) for the oscillator in its steady state, shown 
in Fig.[2]for the model described by Eq. (|2]). It was shown 
in Ref. |6[ that the resonator Wigner functions indicate 
that the crossover is from a fixed-point state in phase 
space to a limit-cycle state. These steady states corre- 
spond classically to a static state, where the amplitude 
is zero, and a state with well-defined amplitude under- 
going harmonic oscillations, respectively. Upon increas- 
ing the coupling strength further, it is possible to drive 
the resonator through a series of first-order crossovers. 
These dynamical crossovers are associated with changes 
between limit cycles of different amplitudes. The first 
order crossovers are also illustrated in Fig. [2j where the 
most-probable resonator state n mp is plotted. This allows 
identification of the mid-points of the crossovers, where 
n mp shows sharp jumps. 




FIG. 2. (color online) A density plot showing the state of 
the resonator in the SSET-resonator system as a function of 
the coupling strength, A, between the resonator and the SSET 
island. Shown is the occupation probability of resonator num- 
ber states P(n) with additional lines showing the expectation 
(n) (solid) and the most probable number state n mp (dashed) . 
The system parameters are T = Vds /eRj, with Rj — h/e 2 , 
Ej = eVds/16 and 7 ex t = 0.0005r. 



B. Quantum Trajectories and the s-Ensemble 

We now turn to the study of the statistics of quantum- 
jump trajectories in the system. A trajectory is a 
time record of quantum jumps associated with particular 
stochastic operators in the master equation [7|, |8|. In the 
SSET-resonator system we have the choice of counting 
charges leaving the SSET island through the right-hand 
lead or counting quanta entering the environment of the 
resonator. The main focus of this paper is on trajectories 
of the resonator, but we study briefly charge trajectories 
in Section [Vj By considering the projection of the den- 
sity matrix p(t) on to the subspace where K events have 
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occurred within the time t, we can define the reduced 
density matrix p^ K \t) such that the probability to ob- 
serve K events in time t is given by Pt(K) = Tr\p( K ^ (t)]. 
After long times, Pt{K) takes a large-deviation form [l4j j: 

P t (K) ~ e - t ^ K/t) (5) 

where cp(K/t) is the large-deviation function which allows 
a complete description of the statistics of K at long times. 
The 8-ensemble 0, H[ enters by considering a moment- 
generating function associated with counting probabili- 
ties 

Zt(s) = Y J Pt(K)e- sK ~e"« (6) 

K 

where the second equality is valid for long times. The 
large deviation functions 6{s) and with k = 

K/t, are related by the Legendre transform 6(s) = 
— mmk[(f(k) + ks] [3, [32|. The s-field is the intensive 
conjugate field to the time-extensive number of events 
K, and the scaled activity k may be used as an order pa- 
rameter to distinguish different dynamical phases. The 
large deviation functions cp(k) and 6(s) jol— Fllj] themselves 
take on roles analogous to those of entropy density and 
free energy density density in equilibrium statistical me- 
chanics. 

The large-deviation function 0(s) will be studied 
throughout the rest of this paper as it encodes the full dis- 
tribution of trajectories Q, [32J. Furthermore, analogous 
to minus a free energy in equilibrium statistical mechan- 
ics, discontinuities in the derivatives of 0(s) have been 
found to correspond to phase transitions in ensembles 
of trajectories of the dynamics Q. The large-deviation 
function 0(s) may be found from a generalized master 
equation [33|, which takes the form of an s-dependent 
modification to Eq. (|2j). First, we introduce an s- 
biased density matrix, p s , defined by 

p s {t) = Y J P (K \t)e- sK . (7) 

K 

Then, when studying the trajectories associated with 
photons leaving the resonator, 6(s) is the largest eigen- 
value of the generalized master equation 

p s = m(Ps) = -\[h,ps] + c™-{ps) + c°r{ps) (8) 

where 

C° S SC (P S ) = Text (e-'ap.a! - ±{<fa,p a }) . (9) 

For s — > 0, the superoperator W s collapses to W and, 
necessarily, 6(s) — > as this corresponds to (the usual 
trace-preserving) physical dynamics. Away from 5 = 0, 
the dynamics are biased by the s-field towards rare tra- 
jectories with, for s > 0, fewer events K within a time 
t or, for s < 0, more events. We will refer to these as, 



respectively, less active and more active rare trajecto- 
ries. While the superoperator W s is no longer trace- 
preserving when 5^0, the mapping is still completely 
positive [34[ as it can straightforwardly be shown to be of 
Kraus form [35|. For all s, effective steady- state proper- 
ties can be determined from the right eigenmatrix of W s 
associated with the largest eigenvalue 0(s). Some peda- 
gogical examples of the s-ensemble applied to quantum 
master equations are discussed in Ref. Q. 



III. NUMERICAL RESULTS 

The SSET-resonator generalized master equation J5J) 
may be expressed in matrix form and diagonalized nu- 
merically, provided the basis for the oscillator is trun- 
cated. For the oscillator damping 7 ex t/r = 0.0005 stud- 
ied in this paper, the Josephson energies will be chosen 
so that the maximum energy of the resonator nfvuj has 
n < 200 [H,|33. However, since the charge on the SSET 
couples to coherences between the oscillator eigenstates, 
it is necessary to keep a basis with off-diagonal density 
matrix elements for the oscillator. In practice, for the 
coupling linear in the oscillator position, it is possible to 
truncate the basis such that only eigenstates which differ 
in energy by mfiw (with m < n) have coherences pre- 
served. The choice of m may be tested numerically to 
ensure that the results are not sensitive to this trunca- 
tion. For the SSET island, coherences are only generated 
between the states |0) and |2), such that propagation of 
the SSET density matrix in time requires inclusion of just 
five of the nine matrix elements. Therefore, the exact 
numerical study of 0(s) requires extraction of the largest 
real eigenvalue of an approximately 5nm x 5nm matrix. 
In the results which follow, we implemented an Arnoldi 
iteration scheme [38] to find 9(s) from the generalized 
master equation. 



A. Resonator Dynamics 

When the SSET bias point is chosen such that AE < 0, 
energy is transferred from the SSET island to the res- 
onator on average. When increasing the strength of the 
coupling between the SSET and the resonator, A, the res- 
onator is driven away from its fixed-point state through a 
continuous crossover into a limit-cycle. Further increases 
in A drive the oscillator through a series of apparently 
first-order dynamical crossovers as the resonator demon- 
strates a series of bistabilities between different limit cy- 
cles, as illustrated in Fig. O 

By measuring the emission of photons from a stripline 
resonator [25|-[28|, it would be possible to infer the state 
of the resonator. Due to the linear coupling to the bath, 
the activity is proportional to both the energy of the 
oscillator hum and the decay rate 7 ex t- The activity may 
be extracted in the s-ensemble from the large deviation 
function 9(s). From Eq. it is clear in the steady state 



that the activity k = (K) /t may be found as: 
1 d 



k = 



t ds 



sK 



K 



s=0 



= -^(o) (10) 



where the prime in the last equality denotes differentia- 
tion with respect to s. Higher derivatives of 0(s) corre- 
spond to higher moments of the distribution of photon 
emissions from the resonator. 

In contrast to studies of full counting statistics j2jj Ho| 
which focus on the s = physical dynamics (see also 
Ref. [41]), we will study the dynamical phases as a func- 
tion of the counting field s. In this way, the dynami- 
cal behavior at s = as, for example, a function of A 
may be understood through its proximity to the phase 
boundaries in the A — s plane. We will first explore the 
full dynamical phase diagram for the case of a resonator 
coupled to an SSET where the quasiparticle decay rate 
is matched to the oscillator frequency, uo = T, and the 
SSET level separation AE is negative, such that there is, 
on average, transfer of energy from the SSET island to 
the resonator. 

For ensembles of trajectories with s ^ 0, we will use 
the activity at non-zero s, found from k(s) = —0 f (s) by 
extending ([TO]) to non-zero s, as the order parameter to 
distinguish dynamical phases. In Fig. [3j k(s) is used 
to construct the A — s phase diagram, shown for SSET 
parameters [3 AE/eV ds = -0.1 and Ej/eV ds = 1/16. 
Figure [3] demonstrates how the signatures of dynamical 
phase transitions in the 8 = physical dynamics shown 
in Fig. [2] may be related to the full A — s phase diagram. 
While at s = we see smooth crossovers in the dynamical 
behavior, for ensembles biased towards rare trajectories 
the crossovers become sharp. The less active trajectories 
for s < show a series of first-order dynamical transitions 
where the activity shows distinct jumps at singular points 
in the A — s plane. For s > 0, the more active trajectories 
undergo a single phase transition provided that s > s c , 
where s c > is the s-coordinate of a critical point at 
which the right-most transition line terminates when A ~ 
0.02. We illustrate the connection between the dynamical 
phase diagram and the 5 = trajectories by plotting 
examples of trajectories at different A in Fig. [3] These 
were found by sampling the full distribution of counting 
events in the non-equilibrium steady state. In particular, 
we note the development of fluctuations in the activity 
on a wide range of timescales close to the critical point 
near A = 0.02. 

The similarity of the dynamical phase diagram in Fig. [3] 
to that derived for the micromaser (Fig. 1 in |8j) further 
demonstrates the parallels between the SSET-resonator 
system and the micromaser noted in [6]. However, the 
SSET-resonator system allows a broader range of dynam- 
ical behavior which can be found by varying its large 
number of tunable parameters. In particular, we can 
vary the strength of coherent oscillations in the SSET by 
altering the Josephson energy Ej, and study the corre- 
sponding dynamical states of the resonator. 
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FIG. 3. (color online) (a) A plot of the activity k as a function 
of coupling strength A and counting field s for system param- 
eters as in Fig. [21 where Ej /eVds = 1/16. With increasing A, 
the less active rare trajectories show a continuous transition 
boundary terminating in a critical point close to s — while 
the more active dynamics exhibit a series of first-order tran- 
sitions which become diffuse as s is increased to zero. Cuts 
through the plot show the activity as a function of A for (b) 
s — —0.1, (c) s = and (d) s — 0.1. Samples of the s = 
quantum-jump trajectories for quanta leaving the resonator 
are shown for couplings (e) A = 0.056, (f) A = 0.041 and for 
(g) A = 0.020 where the system is close to the critical point 
shown in (a). 



In Fig. [H phase diagrams equivalent to Fig. [3] are 
shown for Josephson energies both increased and de- 
creased by a factor of two. For the case of the smaller 
Josephson energy Ej/eV^s = 1/32, the s = dynam- 
ics appears uninteresting. As A is increased from zero, 
the activity shows a weaker version of the continuous 
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FIG. 4. (color online) Plots of the activity k as a function of 
coupling strength A and counting field s for system parameters 
in Fig.[3]but with (a) Ej/eV ds = 1/32 and (b) Ej/eV ds = 1/8 
for fixed F = eVds/h. 



5 = case of the same A. We note that the biased en- 
sembles of trajectories, where 5^0, can be inferred from 
accurate measurements of the distribution of trajectories 
at s = 0: with a knowledge of the large-deviation func- 
tion (p in ([5]), 6(s) can be derived by Legendre transform, 
with k(s) then found by differentiation. In this indirect 
way, the full phase diagram could in principle be gener- 
ated from measurements of the physical system. We are 
therefore able to give a thermodynamic interpretation of 
the distribution of trajectories which gives insights into 
the underlying internal dynamics. However, since the bi- 
ased ensembles may be generated from a knowledge of 
the full distribution of trajectories at s = 0, the fluctu- 
ations at s = must reflect the structure of the phase 
diagram found when favoring rare trajectories. 

We now study the case where Ej is increased from 
the micromaser-like regime illustrated in Fig. [3] We 
show that this introduces new features in the dynamical 
phase diagram even at s = 0. For the Josephson energy 
Ej/eVds = 1/8 shown in Fig.[4j the dynamical transitions 
remain sharp close to s = 0. This is in contrast to Fig. [3] 
where the transition line where s < gradually becomes 
diffuse as it approaches s = 0. In Fig. 2J there is also 
a new crossover in dynamical behavior for s > at cou- 
pling strengths larger than for the initial continuous tran- 
sition. The first-order transition line which approaches 
from s < extends to s = 0, but becomes diffuse for 
positive values of s. These interesting new features go 
beyond what can be seen in the micromaser [8| and their 
origin will be explained using a mean-field treatment in 
Section IIVI In general, we observe that the sharpness of 
the transitions near s — increases with the Josephson 
energy. Correspondingly, the number of quanta in the 
resonator state increases, such that these increases in Ej 
take the system closer to the thermodynamic limit for 
dynamical phase transitions. 



crossover seen for smaller A in Fig. [3^, but none of 
the first-order crossovers. However, the other dynami- 
cal phases which occur when the resonator is oscillating 
in a limit cycle may still be seen when biasing the system 
towards more active rare trajectories. We can interpret 
these findings from the s-ensemble for the physical dy- 
namics as follows: while on average the smaller Joseph- 
son energy restricts the ability of the SSET to populate 
the oscillator with a large number of quanta, on rare oc- 
casions the resonator will be excited to one of the limit 
cycle states where it will release photons into its envi- 
ronment more often. So when biasing the ensemble of 
trajectories increasingly towards those with higher ac- 
tivity, first-order phase boundaries are crossed and we 
may infer that the typical state of the resonator for this 
biased ensemble must be higher in energy than for the 



IV. MEAN-FIELD APPROACH 



In this section we will first present an overview of 
mean- field results [30|, as applied to the SSET- 
resonator system with the Lindblad baths used in this 
paper. These results provide a mean- field theory for the 
resonator amplitude at s = in terms of an effective 
damping, non-linear in the resonator amplitude, arising 
from coupling to the SSET. From this we will develop 
an effective stochastic Liouvillian superoperator for the 
resonator, from which the behavior at s ^ will be de- 
rived. We will then relate these mean-field solutions to 
the exact results found from numerical diagonalization 
of the full generalized master equation (j8j) and explore 
the origin of the dynamical phase boundaries in terms of 
multistabilities in the resonator. 
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A. Mean- field Dynamics at s — 

At s = 0, the dynamical instabilities of the SSET- 
resonator system may be understood from a mean-field 
description of th e dy namics. We follow the approach 
developed in Ref. [30[ . Although the form of the bath for 
the resonator we use differs from that in [30], the final 
mean-field descriptions are the same. 

We define the occupation probabilities of the SSET 
states |0), |1) and |2) respectively by poo, Pn and p 2 2 
and the off-diagonal element of the SSET density matrix 
describing coherences between states |0) and |2) by po2- 
The mean-field equations of motion for the SSET are 
given by 



(Pn) = r((p 2 2> - (pn» 



(11) 



(P22> = -r (p 2 2> + -2^((P02> - (P20» (12) 

(Po2> = (2 (P22} + (pu) - 1) - g (P02) 

+ 1 [A£ + 2a; s C((a) + (a t ))] (p 02 ) (13) 
with the mean-field resonator dynamics described by 



(a) = -iu (a) - ix 8 C((pn)+2 (P22)) 



Text 



(a) (14) 



(ai) = iuj (a^+ix.CUpu) + 2 (p 22 )) - ^fV) .(15) 

The occupation of the SSET state |0) is determined by 
the other two SSET states since (poo) + (pn) + (P22) = 1- 
We have replaced correlators (po2&^) with (P02) ( a ^) 
in ([T3]) . This approximation does not fully capture the 
quantum dynamics, but comparisons with numerics [30| 
show that it provides a reasonably accurate picture of the 
average dynamics of the system, including the limit-cycle 
states of the resonator. 

Because the SSET-resonator coupling is weak, we ex- 
pect the resonator to undergo harmonic oscillations at 
its fundamental frequency uj with an amplitude A which 
changes on time scales much larger than uj~ x . We will 
denote the position about which the resonator oscillates 
(often described as its average or fixed-point position) 
by Xf p , which itself is a function of the SSET param- 
eters and the coupling strength. We will solve for the 
SSET dynamics assuming such harmonic oscillations of 
arbitrary amplitude in the resonator and then use these 
results to find the stable non-equilibrium steady states 
of the resonator. We therefore proceed to solve for the 
SSET dynamics by choosing the following ansatz for the 
resonator state: 



(a) 



raw 
~2h 



(x ip + Ae ut ) 



W) = ] j^(x {p + Ae-« t ). (16) 

Because Ej/hT < 1, we may also consider that the oc- 
cupation probabilities (pn) and (P22) will remain much 



less than unity (30| and can be neglected in (fl3|) . When 
using the ansatz ([T6]h Eq. ([T3]) becomes 



(Pi 



•02/ 



{• 

iEj 
2h 



AE 



2raw x s 
h 



(xf p + A cos ( 



!}(P02> 

(17) 



which can be solved in closed form for a given resonator 
amplitude A. We absorb the A-dependence of (P02) i n 
a time-dependent phase shift a(t) = (2mujAx s /h) smut 
and then find the Fourier series such that (P02) — 
e ia(t) ^ n e iujnt p% 2 [5]. We find Fourier coefficients p% 2 = 
Tp n Jn(— z) , where J n (z) are Bessel functions of the first 
kind, and the parameters ip n are given by 



^ = 



-iEj/h 



2[iun - i(AE/h + 2muj 2 x s Xf p /h - T/2)} ' 



(18) 



We now turn back to the resonator dynamics. With 
our ansatz (fTBjh we use Eqs. (fT4|) and (fT5j) to find an 
equation of motion for the resonator amplitude A. Since 
the amplitude changes on timescales much longer than 
cj _1 , we average A over a time period 2tt/uj to find 



A = 



Text t 1 

■— A ~ UJXs ^ 



2ir/u 



dt((pn) + 2 (p 22 ))sin^ 
(19) 

where the bar on the amplitude signifies that the re- 
lation only holds on time scales long compared with 
The second term on the right-hand side may 
be considered an effective, amplitude-dependent damp- 
ing term Tsset(^4)A but it should be noted that when 
AE < 0, Tsset is actually negative as the SSET pumps 
the resonator. Using Fourier series of the form (pn) = 
E n e ium Vn and (p 22 ) = ]T n e iu;nt p^ 2 we find that 



Tsset (A) A 



((Pii-Pn) + 2(p^ 2 -p 2 - 2 1 )) . 



(20) 

This may be found using the Fourier transforms of 
Eqs. (PU and (H2) and Eq. $1TJ) . Ultimately we find 



ujx s Ej 
Tsset-4 = — Im 



where 



r + iu (r + iu) 2 



(21) 



P = Yi S ^~ mj m+l{z) ~ (^yJm-liz)) J m {z) . 
m 

(22) 

This complicated form of /3 arises from the need to switch 
between the phase-shifted Fourier series for (P02) and the 
Fourier series without phase shifts. The 8 = steady 
states are found as solutions where A = 0. 

As well as finding the steady-state amplitudes of the 
resonator, we may also find Xf p . While the effective 
damping is derived from the lowest frequency oscillat- 
ing terms in the Fourier series for the charge dynamics, 
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the fixed-point position may be found in terms of the 
zero-frequency terms. Using again the Fourier series for 
Eqs. (pi]) and and Eq. ([TTJ), we find 



fffp = -x 8 (Pii + 2^22) = - 



3Ej x s 



2T AE 2 



3E 2 



r 2 



(23) 



In order to find mean-field solutions when s / 0, we 
would need to repeat this analysis with s introduced from 
the start. However, since the full s-dependent master 
equation (j2j) is not trace preserving for general s, the 
method above is not easy to extend to non-zero s. Instead 
we proceed via a different route: we construct an effective 
master equation for the system with a simpler structure. 
This new master equation is chosen to give the mean- 
field results at s = found from the solution above and 
allows solutions with s 7^ to be found using a variational 
method. 



elements of p s to other diagonal elements in the basis of 
number states for the harmonic oscillator, we may con- 
struct from ([24]) a normal operator W s . This operator 
acts on diagonal density matrices in the number basis, 
which we may represent as vectors whose time evolution 
follows: 

p s = W s p s = |e" s 7extV / ata + la - Text^a 

+ g(n)\/a)a + la 1 " - g(n){a) a + 1) jp fl . (27) 

Using a variational ansatz of coherent states as in [gj], 
we set a = e z6 \fn and a) = e~ %6 \fn. We can then find 
a mean-field estimate of the largest eigenvalue, 0(s), by 
maximizing W s with respect to 5 and n. This is done by 
finding the values of S and n which satisfy dW s /dS = 
and dW s /dn = 0. The first equation allows S to be 
eliminated so that 



B. Effective Stochastic Master Equation 

Using the above result for the effective (negative) 
damping term for the resonator, arising from the weak 
coupling to the SSET island, we now construct an effec- 
tive master equation for the resonator. By construction, 
this master equation will correctly describe the average 
s = dynamics to within the mean-field approximations 
above. However, noise arising from the stochastic quasi- 
particle decay process cannot be captured. The structure 
of the master equation is simpler as it just involves two 
Lindblad operators; one is associated with the environ- 
ment of the resonator and the other describes the driving 
due to the SSET. The complexity of the mean field so- 
lution will be contained in these driving terms, which 
must be amplitude dependent. By applying the counting 
field s to the trajectories defined by the time record of 
quanta entering the environment from the oscillator, this 
generalized master equation is 

p s = Ws(ps) = Lpstf - ^{L^L.ps} 

+ Text (e~ s ap s a ] - ^{a f a, p s }^ (24) 
where the new Lindblad operator is defined by 

L=^g{n)a ] (25) 

where n = a) a and 



s/2 A 



I n g(n) 

Text 



(28) 



?(n) = -7ssetW ■ 



(26) 



The amplitude-dependent negative damping allows con- 
struction of an n-dependent driving term, where n is 
the number of quanta in the oscillator state, which is 
related to its amplitude via n = A 2 (muj/2h). In addi- 
tion to its dependence on n, g(n) depends implicitly on 
the SSET parameters and the coupling strength, A. Be- 
cause this stochastic description only couples diagonal 



and a) = n/a. Inserting these into W s yields 



W s =2e 3/2 y/n(n + 1)^9 {n) - Text^ - g{n){n + 1) 

(29) 

which may be maximized with respect to n numerically. 

For the SSET parameters which show the micromaser- 
like behavior demonstrated in Fig. [3j the mean-field solu- 
tion found from Eq. ([29]) is shown in Fig. [5] The case of 
larger Ej /eVds = 1/8 shown in Fig. [4] may be compared 
with the mean-field results in Fig. [6] For larger Ej 
there is excellent agreement, but for lower Ej values the 
phase diagrams, though qualitatively similar, differ quite 
noticeably in the exact locations of transitions. (Accord- 
ingly, the mean-field result captures the dynamical phase 
diagram less accurately when Ej/eV^s — 1/32.) We sug- 
gest that this is most likely due to the SSET charge 
fluctuations which are neglected in Eq. ([24]) . In the pa- 
rameter regime studied here, these fluctuations are more 
significant when Ej is smaller. Support for this view 
comes from the numerical result for the full distribution 
P(n) for the occupation of the resonator number states 
for Ej/eV ds = 1/8 shown in Fig. [7(b). This should 
be contrasted with the P(n) distribution in Fig. [2] for 
Ej/eVds — 1/16, which shows significant blurring be- 
tween the peaks in n corresponding to different limit cy- 
cles, when compared to Fig. [7(b). These effects will not 
be captured by the mean-field solution as it neglects these 
fluctuations. 

Another observation we make is that the mean-field 
phase diagrams show an extra phase boundary when 
s > for Ej/eVds = 1/8, which is also present in the 
numerical result in Fig. [4j We will show that the appear- 
ance of this new phase boundary may be understood in 
terms of the mean-field equations involving the non-linear 
driving term due to the SSET. 
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C. Phase Transitions at s — 



At s = 0, the dynamical states of the resonator are a 
set of limit cycles whose amplitudes are determined by 
balancing the effects of the dissipation due to the envi- 
ronment, which is linear in the resonator amplitude, and 
the non-linear driving due to the SSET. This is clear from 
Eq. ([29]) when written in the form 



0.12 



(V7ext^ - 

+2(e- s / 2 



V#(n)(n + 1)) 2 



l)Vrc(rc + l)7ext0(n). (30) 



At s = we require 6(0) = 0. This occurs when Eq. ([30]) 
is maximized since, for s = 0, W s < 0. Multist abilities 
for the resonator state occur when there is more than one 
value of n which satisfies W s = 0. From Eq. ([30]) . we can 
also see that these multistabilities are associated with 
phase boundaries at s = 0. To see this, consider the case 
where many values n^, with n^+i > n^, satisfy W s = 
with s = 0. Now, noting that g(n) is a smooth function 
of n, if we perturb s from zero by a small amount 5s, 
changes in the positions of maxima rti will be of order Ss 
and, therefore, much smaller than n^ + i — rti. Now if we 
consider the values of W s when introducing small shifts 
Ss at the maxima we see that 



W Ss = -^ v / 7 ext ^(n i )n i (n i + 1) = Ssj e: 



(31) 



to first order in Ss. From Eq. ([3T]) we can see that if 
Ss > 0, maximization of W s is found by selecting the 
smallest n^. Conversely, for Ss < 0, W s is maximized by 
picking the largest rti. Therefore we understand multi- 
stabilities, where there are two or more solutions to the 




FIG. 5. (color online) Mean field: A plot of the activity k as a 
function of coupling strength A and counting field s for system 
parameters 7 ex t/r = 0.0005, u/F = 1, AE/eV ds = -0.1, 
Ej/eV ds = 1/16 for fixed F = eV ds /h. 
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FIG. 6. (color online) Mean field: A plot of the activity k as a 
function of coupling strength A and counting field s for system 
parameters 7 ex t/r = 0.0005, uj/F = 1, AE/eV ds = -0.1, 
Ej/eV ds = 1/8 for fixed F = eV ds /h. 



s = mean-field equations, are associated with s = 
mean-field dynamical phase boundaries. Upon changing 
A, phase boundaries extend into the s > or s < half- 
planes when, respectively, the smallest rti solution or the 
largest rti solution changes. With this insight, we now 
turn to the question of why increasing Ej introduces a 
new phase boundary on the less active side of the phase 
diagram. 

For Ej/eVds = 1/16, as shown in Fig. [2j the onset of 
bistability occurs at A ~ 0.012 and the resonator remains 
multist able as A is increased further. In contrast, when 
Ej/eVds = 1/8? there exists a second region of monos- 
tablilty as A is increased beyond the bistable region. In 
Fig.[7J we show the full P(n) distribution for the occupa- 
tion of the oscillator number states when Ej/eVa s = 1/8? 
obtained from exact diagonalization of the generalized 
master operator. The initial onset of bistability occurs 
at A ^ 0.06 with monostability again beyond A ~ 0.077. 
This is consistent with the phase boundary shown at 
s = between these values of A in both the exact phase 
diagram in Fig. 2] and the mean-field phase diagram in 
Fig. [6l In Fig. [7J we also show that for s perturbed 
slightly from the zero, the bistability disappears from 
the full P(n) distribution as one of the two possible limit 
cycles is selected, consistent with Eq. ([30]) . 

We now turn specifically to the origin of the second 
monostability present for Ej/eVds = 1/8. We may un- 
derstand this with reference to the effective negative 
damping 7sset on the oscillator provided by the SSET. 
If we neglect the effects of the SSET on the fixed-point 
displacement of the resonator, Xf p , from Eqs. ([T8]h ([2T]) 
and ([22]) . we note that 7sset^4/^s scales with Ej so 
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jext/EjX 2 . To understand the existence of new limit 
cycles upon increasing A, it should be noted that these 
occur when A is large enough that such lines have shallow 
enough slope to intersect peaks in the effective damping 
at larger X 2 A/x s . The case where j ext /EjX 2 = 1.8 is 
plotted to illustrate that while monostability exists be- 
yond the first region of bistability when Ej/eVd s = 1/8, 
for Ej/eVds = 1/16 these solutions are not possible due 
to the smaller ratio Xf p /x s . 
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FIG. 8. Plots of the effective negative damping on the 
resonator due to the charge on the SSET island. Plotted 
are the damping curves corresponding to Josephson energies 
Ej/eVds = 1/16 (solid) and Ej/eV ds = 1/8 (dashed), with a 
magnified plot of the area of interest shown inset. All other 
parameters are as in Figs. [5] and [6] The dotted line of gradi- 
ent 7ext/^jA 2 = 1.8 illustrates the origin of the second region 
of monostability when Ej/eVds = 1/8 since this line has just 
one intersection with the dashed curve. 



FIG. 7. (color online) The state of the resonator close to 
s = 0. Full distributions of P(n) are plotted as a function of 
A for the parameters in Fig. [4] with Ej/VdS = 1/8. Shown are 
three s values: (a) s = —0.005, (b) s = and (c) s = 0.005. 
Shown with solid lines are the expectation values (n). 



that 7sset^4/^j^s is a function of XA/x s . Therefore, 
where Xf p may be neglected, we expect the same struc- 
ture of limit cycles and multistabilities at different Ej 
from the mean- field theory. However, if we include the 
effects of the modified fixed point using Eq. ([23]) , we find 
that the effective damping is modified by a small but 
significant amount. In Fig. [8j we show how Xf p per- 
turbs the form of 7sset upon increasing Ej. By plot- 
ting — 7sset^4/^j^s against XA/x Sl we show disappear- 
ance of bistability when increasing Ej/eVd s from 1/16 to 
1/8. Limit cycle solutions occur where — 7sset = 7ext- 
These solutions are illustrated on Fig. [8] by the intersec- 
tion of — jsset A/ EjX s with lines of constant gradient 



V. QUASIPARTICLE TRAJECTORIES 

We now turn to the trajectories of quasiparticles de- 
caying into the right-hand lead of the SSET. For oscilla- 
tors with uj/T < 1, there are no features in the current 
corresponding to changes in the dynamical state of the 
resonator. So far we have considered the case of a res- 
onator with relative frequency uj/T = 1. In this regime, 
we are not able to resolve individual features in the res- 
onator energy as a function of AE, which arise from in- 
elastic absorption of energy quanta from the SSET when 
AE = —jhw, for integral j > 0. However, if we con- 
sider a fast oscillator where uj/T ^> 1, the timescales for 
the SSET and resonator dynamics become separated and 
the oscillator does show distinct peaks in energy when 
AE = —jhuj [6]. Between these peaks, the SSET has lit- 
tle influence on the resonator dynamics. In Fig.[9j we plot 
numerical results for the AE — s phase diagram in this 
regime. The resonator absorption peaks separated by 
Huj correspond to strong features in the phase diagram. 
We have also applied the mean field theory developed in 
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Section IIVI to this parameter regime. This method also 
exhibits features separated by two but only shows quali- 
tative agreement with Fig. [9j a result which we attribute 
to the small resonator energy across much of the AE — s 
phase diagram in this case. 




FIG. 9. (color online) A plot of the AE — s phase dia- 
gram showing the activity for quantum jumps associated with 
the environment of the resonator. The vertical axis shows 
AE/hw. The system parameters are uj/T — 3.5, 7ext/r = 
0.002, A = 0.2 with Josephson energy Ej/V ds = 1/16. 



The quasiparticle current in the resolved side-band 
limit, occurring when uj/T ^> 1, does exhibit maxima 
which coincide with peaks in the energy of the resonator 
as a function of AE [36[. We understand the emer- 
gence of these peaks as follows: when uj/T ^> 1, the 
separation of timescales of resonator and SSET means 
that for most bias points, the quasiparticle current is low 
and not strongly influenced by the presence of the res- 
onator. However, at the particular values AE = —jhuj, 
absorption of energy by the resonator from the charges 
occurs allowing a greatly enhanced charge flow through 
the SSET island (analogous to an inelastic tunneling cur- 
rent). This correlation between resonator energy and 
charge dynamics motivates our study of the trajectories 
of SSET quasiparticles. 

We now apply the s-ensemble to charges leaving the 
SSET island. We wish to explore the extent to which the 
trajectories of the SSET quasiparticles mirror the tra- 
jectories of the photons entering the environment of the 
resonator. To do this, we apply the s-ensemble to the 
counting of charges leaving the SSET island by diagonal- 



izing the superoperator W^ p (p) defined by 

(p) = -{ i H > P] + C P - (P) + £° sc - (P) (32) 

where 

CriP) =re- s (|l)(2|p|2)(l| + |0)(l|p|l)(0|) 

-£({|2)<2|,p} + {|l)(l|,p}). (33) 

The large-deviation function for quasiparticle trajecto- 
ries, #q. p .(s), is given by the largest real eigenvalue of 
W^ p (p). We will study the dynamics of SSET charges 
using the quasiparticle activity — d# q . p .(s)/ds. 

For a resonator where uj/T < 1, we see no transitions in 
the dynamical behavior of the SSET charges, consistent 
with the known s = behavior [36|. However, for the 
system parameters in Fig. [9] corresponding to a fast os- 
cillator, we find that the AE — s phase diagram, shown in 
Fig. [TUJ shows similar structure to the dynamical phase 
diagram for the trajectories of quanta leaving the res- 
onator in Fig. [9] In the resolved side bands, where mul- 
tiples of the oscillator level spacing closely match the level 
spacing of the SSET charge states, the inelastic interac- 
tion between SSET and resonator enhances the rate at 
which quasiparticles decay. Therefore the trajectories of 
both resonator photons and charge quasiparticles allow 
inference of the resonator state in this regime. However, 
the peaks in the respective activities in Figs. l9land[T0lat 
AE = show a marked difference. When AE = 0, there 
is no significant driving of the resonator by the SSET and 
so the number of quanta leaving the resonator is small. 
In contrast, the decay of quasiparticles into the right- 
hand lead is large, such that the activity when counting 
charge quanta is large at AE = 0. Finally, we mention 
another interesting feature of Fig. [9] When biasing to- 
wards more-active trajectories, we find that the activity 
associated with counting quanta from the resonator is 
bigger when AE = —2huj than when AE = —huj. This 
effect is only seen for biased ensembles with s < 0. 



VI. CONCLUSIONS 

In this paper we have used the s-ensemble approach to 
study the dynamics of a complex system with two com- 
ponents, each of which has non-trivial dynamics. The 
coupling of an SSET and a resonator has been studied 
using mean-field techniques previously [30j |. We devised 
a method to extend a standard mean-field theory for 
the coupled system to enable us to study ensembles of 
trajectories with non-zero s. Comparing the results of 
our mean-field theory to those employing essentially ex- 
act numerical diagonalization of the generalized master 
equation, we found good agreement for system parame- 
ters where the resonator is driven into high energy states 
of self-sustaining oscillation. 

We find the accuracy of mean-field results improves 
for large Josephson energy, where the resonator energy 
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FIG. 10. (color online) A plot of the AE — s phase diagram 
using the quasiparticle activity, found from the generalised 
master equation (|32p . The vertical axis shows AE/ftw. The 
system parameters are identical to those in Fig. [9] 



is large. We were able to use our method to understand 
the correspondence between multistabilities in the state 
of the resonator and boundaries in the dynamical phase 
diagram at s = 0. This demonstrates how the dynamical 
phase diagram for trajectories can be used to infer prop- 
erties of the internal dynamics of the system. However, 
even when there are no phase boundaries near 5 = 0, 
as we showed to be the case for small Josephson energy, 
signatures of the underlying complex internal dynamics 
are encoded in rare trajectories, which in the s-ensemble 
correspond to dynamical phases away from 5 = 0. Thus 
we demonstrated how the s-ensemble provides a method 
for interpretation of accurate measurements of the distri- 
bution of trajectories. 

Exploring the SSET-resonator system also allowed us 
to examine trajectories in the s-ensemble created from 
different operators. As each of the coupled components 
is itself an open system, we were able to study both tra- 
jectories formed by measuring charges in the right-hand 
SSET lead and photons emitted by the resonator. Pre- 
vious studies have shown that fast oscillators imprint a 
signature of their energy absorption from the SSET on 
the SSET charge statistics [36|. We demonstrated that, 
in this regime, the phase structure of trajectories of de- 
caying quasiparticles mirrors that found by measuring 
photons emitted by the resonator. 

Finally, we emphasize the potential significance of 
this work with regard to future experiments. SSET- 
resonator systems have been constructed with a stripline 



resonator [24|. Impressive progr ess towards measuring 
single microwave photons |25l-[28| makes measuring the 
trajectories studied in this work an exciting possibility. 
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Appendix: Practical Requirements 

In this Appendix we give a brief overview of the practi- 
cal requirements for reaching a regime where the dynam- 
ics of a SSET-resonator system is described by Eqs. (pQ) 
to (jlj). Fuller accounts are available in works which fo- 
cus on the theoryof SSETs [H, or SSET-resonator 
systems [H, Hfll Hll |42| and in experimental studies of 
SSETs coupled to either a mechanical [23| or electrical 
resonator [241 ]. 

We start by considering the conditions that must be 
fulfilled for the current through the SSET to be domi- 
nated by JQP processes, which are composed of cycles 
involving a combination of coherent Cooper-pair tunnel- 
ing and quasiparticle decay. The JQP resonances occur 
in the low temperature regime where the electrostatic 
charging energy of the SSET island, E c = e 2 /2Cs, with 
the total capacitance of the island, is much larger 
than k^T. Within this regime, the number of differ- 
ent charge states accessible to the SSET island becomes 
severely limited and we can treat it as a few-level system. 
The other two important energy scales for the SSET are 
the superconducting gap, A, and the Josephson energy 
associated with the junctions Ej = hA/(8e 2 Rj), where 
Rj is the resistance of the Josephson junctions. 

A JQP resonance occurs when the gate and bias volt- 
ages, V g and Vds respectively, are chosen so that the elec- 
trostatic energy difference between states on the island 
differing by one Cooper-pair, AE, is zero and the drain- 
source voltage is large enough so that both of the quasi- 
particle tunneling processes associated with the cycle are 
allowed. The gate voltage induces a polarization charge 
n g = (CgVg + CjVds)/e on the island [HI, where C g and 
Cj are the capacitances of the gate and the Josephson 
junctions respectively. This then leads to an electrostatic 
energy difference AE = —4E c (n g — n — 1) where n is an 
integer, corresponding to the number of Cooper-pairs on 
the island, (we take n = for simplicity in the main text) 
and hence resonance occurs along lines in the V g — Vds 
plane given by n g = n + 1 [18[- Both of the quasipar- 
ticle processes involved in the JQP cycle are allowed for 
eVds > 2 A + E c ; provided the charging energy is large 
enough compared to the gap, E c > 2 A/3, the JQP cy- 
cle dominates the current up to eVds = 4 A, at which 
point current can flow through the motion of quasiparti- 
cles alone. The master equation description of the quasi- 
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particle tunneling processes is derived assuming high re- 
sistances for the junctions [18l-l20|, Rj ^> h/e 2 . 

The coupling of a nanomechanical resonator to a SSET 
is discussed in detail in [18|, [20| , while the case of cou- 
pling to a superconducting stripline resonator is dis- 
cussed in [24[, though from a theoretical point of view 
the form of the coupling is essentially the same in both 
cases. In the case of a mechanical resonator, a metal 
layer is added to a nanomechanical beam fabricated from 
a semiconductor which is adjacent to the SSET island. 
The beam acts as an additional voltage gate and the 



length-scale, x s , which describes the coupling between 
the island and the resonator [see Eq. JTJ] is given by 
x s = 2E c C rn V rn /(emuj 2 d), where C m is the beam-island 
capacitance, V m is the voltage applied to the beam, d 
is the distance between the beam and the island and 
m the effective mass of the resonator. Temperatures 
~ 30 mK [23j are routinely used in experiments with 
SSETs. We note that experiments using a superconduct- 
ing stripline resonator have the advantage that the rele- 
vant mode frequency of the stripline, uj are large enough 
so that fux) ^> k^T under typical conditions. For example, 
in the experiment using a stripline resonator reported in 
Ref. [24] the mode frequency was u/2tt = 9.9 GHz. 
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